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Macroscopic properties of inultiferroics, the systems that show si- 
multaneously two types of ordering, could be controlled by the 
external fields of different nature. We analyze the behavior of mul- 
tiferroics with antiferro-(AFM) and ferromagnetic (FM) ordering 
under the action of external magnetic and stress fields. A combi- 
nation of these two fields makes it possible to achieve macroscopic 
states with different domain structures. The two-domain state ob- 
tained in this way shows a linear dependence of macroscopic strain 
vs magnetic field which is unusual for AFMs. A small but nonzero 
stress applied to the sample can also result in the bias of the mag- 
netization vs magnetic field dependence. 



1. Introduction 

Present-day technologies of information storage are 
based on the i) combination of the materials with differ- 
ent magnetic, electronic, elastic properties; ii) step-like 
variation of the magnetic, electronic, elastic properties 
within a small, nano-sized sample; in) combined appli- 
cation of different fields to implement the control over 
macroscopic properties. The first tendency initiated the 
growing interest in the natural and synthetic multifer- 
roics, i.e., materials that show simultaneously at least 
two different order parameters (like the magnetic and 
electric polarization) . The second tendency is related to 
the appearance of a giant variation (as compared with 
that in homogeneous samples) of macroscopic proper- 
ties like the magnetoresistance and, as a consequence, 
the high susceptibility of samples to external fields. The 



last approach opens a way to produce macroscopic states 
unreacheable in other ways or the desired states much 
more effectively. It should be noted that the sensitivity 
of multiferroics to two fields of different nature instead 
of one should be much greater than that of the single- 
ferroics^. Really, in multiferroics, two order parameters 
have different origin and different symmetry and even 
may appear at different temperatures, while the second- 
order parameter in single-ferroics, which is also the sec- 
ondary one, contributes to the values of characteristic 
fields and susceptibilities, but gives rise to no new qual- 
itative effects. So, we can expect that a combination of 
two fields in multiferroics results in the variety of states 
with different symmetries that could be used for the ef- 
fective information storage. 

Many multiferroics show strong magnetoelastic effects 
along with ferro- or antiferromagnetic (FM and AFM, 
respectively) and ferroelectric ordering (see, e.g., [1-5]). 
In this case, the macroscopic properties could be con- 
trolled by the simultaneous application of stress and 
magnetic/electric fields. An external stress can be used 
for the compensation of internal strains that produce 
the undesirable backswitching of a ferroelectric polariza- 
tion in some multiferroics (like BiFeOa) [4J. Additional 
stresses could be also produced in the course of epitaxial 
growth (on a substrate with properly chosen misfit), like 
it was done, e.g., in Ref. [6]. Such a "plugged" stress 



^ Single-ferroics could also be controlled by two different fields in 
the case where they show some cross-correlated effects like the 
magnetostriction or the piezoeffect. 
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can stabilize a certain type of domain structure (e.g., 
by eliminating undesirable domains) and thus facilitate 
the control over macroscopic properties by external mag- 
netic/electric fields. 

A combination of stress and magnetic/electric fields 
is also interesting from the fundamental point of view, 
because these fields have not only different symmetries, 
but also different tensor nature. Here, we study the be- 
havior of a multidomain magnetic system in the pres- 
ence of stress and magnetic fields. A good example of 
such a system can be found among materials that show 
simultaneously the FM and AFM orderings in various 
systems of magnetic ions (see, e.g., Mn-based antiper- 
ovskites [7, 8] and a high-temperature superconductor 
Sr2Cu304Cl2 [9]) or in the different parts of a sample 
(like manganite Pro.sCao.sMnOa [2, 3], in which the FM 
and AFM phases coexist in a wide range of the temper- 
ature and the field). It is also important that an AFM 
subsystem is usually a source of the pronounced magne- 
tostrictive strain and thus interacts with a stress, while 
the FM subsystem is strongly coupled with an exter- 
nal magnetic field. Moreover, in the finite-size samples, 
the AFM-induced magnetostriction and the FM magne- 
tization are the sources of different shape-induced (de- 
stressing and demagnetizing) effects that compete with 
external fields and may also compete with each other. 

The main question which we address below is whether 
a combination of magnetic and stress fields may reveal 
some new features and facilitate the control over macro- 
scopic properties of a sample. We will analyze the phase 
diagrams of the homogeneous and multidomain states 
and the field dependence of the macroscopic strain and 
the magnetization. We argue that the presence of the 
FM component gives rise to a substantial increase in the 
strain susceptibility to the external magnetic field, and 
the application of a stress can induce the magnetic bias 
seen as a shift of the magnetization curve. 

2. Combined Influence of a Magnetic Field and 
a Stress on the Magnetic System: Qualitative 
Considerations 

Let us consider a system that shows simultaneously FM 
and AFM orderings, no matter how it is realized - as a 
mixture of phases (like in perovskite manganites [2, 3]) or 
as a FM+AFM multiferroic (like Sr2Cu304Cl2 discussed 
in details below) . Suppose further that the system shows 
a noticeable magnetoelastic effect which originates from 
the coupling between the AFM order parameter and the 
lattice. Which mechanisms are responsible for the equi- 
librium domain structure (DS) in this case, and how this 



DS could be controlled? It is well known that the frac- 
tion of FM domains (in a pure FM [10] or in a mixture 
with the AFM phase [II, I2[) is governed by magne- 
tostatic effects usually described with the stray energy 
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where M™^'^''° is the macroscopic magnetization, brack- 
ets (...) mean the averaging over the sample volume 
y, and the second-rank demagnetization tensor N'^™ de- 
pends upon the sample shape. 

On the contrary, the AFM vectors themselves do not 
produce any magnetostatic field. However, due to mag- 
netostriction, an AFM component produces the so-called 
destressing fields of elastic nature [I3[ that give rise to 
the formation of an equilibrium DS. 

In the case of an FM+AFM phase mixture, the in- 
clusions of a new (FM or AFM, depending upon the 
direction of a phase transition) phase could be consid- 
ered as elastic dipoles (due to the strain misfit between 
the FM and AFM phases). In the case of AFM+FM 
multiferroic, the analogous dipoles arise at the sample 
surface (see [14] for details). In both cases, the elas- 
tic dipoles produce long-range fields that could be taken 
into account phenomenologically with the use of the de- 
stressing energy 
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where L is an AFM order parameter, and the 4-th rank 
destressing tensor Nj^i^, hke A(^™, depends upon the 
shape of the sample or/and inclusions and upon the mag- 
netoelastic coupling strength. 

Both the stray (Eq. (1)) and destressing (Eq. (2)) 
contributions to the free energy of a sample constrain a 
variation of the domain fraction in the presence of ex- 
ternal fields. Though both contributions originate from 
different physical mechanisms and even have different 
symmetry tensor characters, their values could be com- 
parable in general case and should be taken into account 
on the same footing. 

The infiuences of a magnetic field and a stress on the 
DS of the FM+AFM system are also different. The ex- 
ternal magnetic field H discerns the domains with oppo- 
site orientations of the macroscopic magnetization Mp 
of FM (or FM component) and also the domains with 
the perpendicular orientation of the AFM vector L, as 
follows from the standard expression for the Zeeman en- 
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ergy^ (per unit volume): 



WZo 



X, 



-(M'"''"'°)H--([L xH]^), 



(3) 



where x is the magnetic susceptibility of the AFM sub- 
system, and the macroscopic magnetization M'"^'^''° in- 
cludes a contribution from both FM and AFM subsys- 
tems, see below, Eq. (6). It should be stressed that the 
domains with L||H are always unfavorable in a magnetic 
field and move away under the field action. 

On the contrary, the mechanical stress (described by 
the symmetric tensor ajt) has no effect on the FM com- 
ponent (unless it produces the magnetostriction) but 
may discern the AFM domains with the perpendicular 
orientation of L, as can be deduced from the expression 
for the magnetoelastic energy (where the elastic isotropy 
is supposed): 
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where A is the magnetoelastic constant, Cjk are the elas- 
tic moduli, and we have omitted the terms related to 
the isotropic (hydrostatic) pressure as immaterial for the 
further consideration. 

If, for example, a tensile/compressive stress is applied 
along the n||H axis, then the combined action of the 
magnetic field and the stress is described by the term 
(obtained from (3) and (4)): 



W'ficld 



-(M'°^"°)H 



1 



xH' + —){[Lxn]'), (5) 
C44/ 



Tr(7/3 



where we introduced the value of stress a 
nan. 

It is quite obvious from Eq. (5) that the action of a 
magnetic field on the AFM subsystem is equivalent to 
the action of a stress with fixed sign. So, a combination 
of H and a gives possibility to enhance (if Act > 0) or 
wipe out (if Act oc —C44xH^) the inequivalence of differ- 
ent AFM domains. Thus, it follows from general consid- 
erations that the combined application of a stress and a 
magnetic field opens a way to act selectively on that or 
those types of domains and to control macroscopic prop- 
erties of the sample (such as the magnetization, elonga- 
tion, magnetoresistance, etc.). 

In the next sections, we consider the typical effects 
produced by two fields in a relatively simple but non- 
trivial AFM+FM muhiferroic Sr2Cu304Cl2. 



■^ We distinguish a magnetization of the FM component Mp from 
a noncompensated magnetization M of AFM. The later results 
from the field-induced canting of the magnetic sublattices of 
AFM. 




X, X 



y y 

Domain A Domain B 
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Fig. 1. (Color online) Magnetic structure of CU3O4 layer in two 
different configurations (domains). Magnetic field is parallel to 
(110). Two types of magnetic ions are represented with filled and 
hollow circles. FM ordered moments of Cull could be (a) parallel 
(domain A) or (6) perpendicular (domain B) to the applied mag- 
netic field. Small canting of the Cul spins induced by the external 
magnetic field is ignored 



3. Model 

The crystal structure of high-temperature superconduct- 
ing cuprates Sr2Cu304Cl2 and Ba2Cu304Cl2 consists of 
CU3O4 planes separated by spacer layers of SrCl or BaCl 
[15, 17]. Two types of magnetic ions, Cul and Cull (see 
Fig. 1) form two interpenetrating square lattices within 
CU3O4 planes. 

Within the temperature interval Tu w 40 K< T < 
Ti sa 380 K, the ions of the first type (Cul) are AFM 
ordered, while the ions of the second type (Cull) bear a 
small but nonzero FM moment'^. According to the ex- 
periments [18], the mutual orientation of Cul and Cull 
moments depends upon the direction of the external 
magnetic field and can be either perpendicular or par- 
allel. In other words, the magnetic structure consists 
of two weakly coupled subsystems, namely, the AFM 
subsystem localized on Cul ions and the FM one local- 
ized on Cull ions. The FM subsystem is unambigu- 
ously described by the magnetization vector Mp, and 
the AFM subsystem is described by two vectors: AFM 
vector L = (Si — S2 + S3 — 84)74 and ferromagnetic 
vector M = ^ Sj/4 (numeration of Cul sites is shown 
in Fig. 1). The macroscopic magnetization M'"'^'^''" (per 
unit volume) is defined as a sum of FM and AFM mag- 



^ According to Ref. [18], the FM moments at Cull ions result 
from the anisotropic "pseudodipolar" interactions between Cul 
and Cull. 
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Fig. 2. (Color online) Four types of magnetic domains. Axes x and 
y are parallel to (100) crystal directions. The external magnetic 
field H||[110] (if any). Types A and B have different orientations of 
the AFM vector, types 1 and 2 correspond to opposite directions 
of the FM vector Mp. Ellipse (dash line) images the supposed 
shape of the sample and its orientation (axes X, Y) with respect 
to crystal axes 



netizations as follows: 



M" 



= Mf 



M. 



(6) 



In the absence of an external field, the FM moments at 
Cull sites are oriented along (f 10) crystal directions per- 
pendicular to the staggered magnetizations of the AFM 
subsystem, as shown in Fig. I. Due to the tetragonal 
symmetry of the crystal (space group lA/mmm), an 
equilibrium magnetic structure can be realized in four 
types of equivalent domains, as shown in Figs. 1 and 
2. Domains of type A and B could be thought of as 
AFM domains, because they correspond to different ori- 
entations of the L vector and thus are sensitive to the 
orientation of the magnetic field H or the stress a with 
respect to the crystal axes (see Fig. 2). Types Al and 
A2 (and, respectively, Bl and B2) are FM domains; they 
have opposite directions of the Mp vector and could be 
removed from the sample by H||Mf. 

If the external (magnetic and stress) fields are ap- 
plied within the xy (and, equivalently, XY) plane and 
their values are small enough as compared with the ex- 
change interactions between different Cu sites, the mag- 
netic structure of Sr2Cu304Cl2 can be unambiguously 



described with two angle variables, as shown in Fig. 2: 
Lx — Mq cos 9, Ly ~ Mq sin 0] 



My 



m-pMo cos if, Mpy — TOpA/o sin ip. 



(7) 



Here, mF(=10^'^ for Sr2Cu304Cl2 [18]) is a dimension- 
less constant that represents the ratio between the spin 
moments localized on Cull and Cul sites. We have also 
taken into account that, far below the Neel temperature, 
the values of sublattice magnetizations Mq = 27.4 Gs 
(that corresponds to the spin s = 1/2 per Cul site) and 
Mp are saturated and constant. 

The phenomenological description of the DS is based 
on the analysis of the free energy potential <& of the 
sample. We consider four constituents of $: mag- 
netic $mag, with account of magnetoelastic interactions; 
shape-dependent stray (demagnetizing) Eq. (1) and de- 
stressing Eq. (2) energies, and the energy $fioid = 
^Wficid (see Eq. (5)) of the external fields, as explained 
above: 
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The magnetic energy of Sr2Cu304Cl2 crystal in the 
mean field approximation is well established [9, 16, 18] 
and was analyzed in details, along with $stray and $dostj 
in ]19]. Here, we start directly from the simplified ex- 
pression for the specific potential <j) = ^/V expressed 
through the angles 9 and ip: 



(j) = 4JpdmF(cos(6' + ip)) + K\\ { 
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+AN'"'^sin2{e-^)y 



(9) 



In (9), ip is an angle between the unit vector nl|H 
and the a;-axis, a mechanical stress is applied along n as 
described above, N^'^ are the components of the demag- 
netization, and N'^'^^, AN'^°^ are those of the destressing 
tensor which are important for the chosen geometry of 
the sample (thin pillar of the thickness c with an elliptic 
cross-section, whose principal axes X and Y are parallel 
to (110) directions within the CU3O4 layers, see Fig. 2), 
and N2°^ is a shape-induced anisotropy of the AFM sub- 
system. The meaning and values (in Oe) of phenomeno- 
logical constants used in (9) are given in Table 1. Here 
and for the rest of the paper, we use the values in Oe 
instead of energy units (say, (j) — > (J>/Mq, etc.). For the 
rest of the paper, we assume that the magnetic field and 
the stress are applied in parallel to one of the principal 
axes of the sample that coincides with an easy axis for 
the AFM vector, namely. n\\X, (see Fig. 2), so ^/i = 7r/4. 
This situation corresponds to the experimental setup in 
Refs. [9, 16]. 

4. Homogeneous States in the Presence of Two 
Fields 

Let us start from the analysis of possible homogeneous 
states that could be realized in the thermodynamic limit 
(infinite sample, all TVdm^^dos ^ g^^ i^ ^.j^jg ^^^^^^ ^^^ 

minimization of (j) with respect to the magnetic variables 
9 and ip gives rise to four solutions labeled as Al,2 and 
Bl,2 (see Fig. 2). Corresponding equilibrium values at 
H = {] and cr = are 



Al: 0A1 



-tt/A, ipAi = 7r/4; 



Bl: 0B1 = 5^/4, ^Bi = 37r/4; 
A2: 6a2 = 37r/4, ^a2 = 57r/4; 

B2: 6Ib2 = 7r/4, v3b2 = -7r/4. 



(10) 



Table 1. Parameters used in the free energy potential 
[Eq. (9)]. Raw data in meV taken from [9, 17, 18] are 
converted to Oe 



Parameter 



Meaning 



Value (Oe) 



Jo 

Jav 
Jpd 

K„ 



Cul-Cul superexchange (in-plane) 1.02 X 10^ 

isotropic pseudodipolar interaction —9.4 X 10^ 

anisotropic pseudodipolar interaction —2.1 X 10^ 




Fig. 3. (Color online) Stability ranges of different homogeneous 
states in the presence of a field H and a stress a. Solid lines bound 
the stability ranges of A2 (line 1), B (lines 2, 3) and Al (line 4) 
states 

It should be stressed that, in contrast to pure AFMs, 
the configurations with (Mf,L) and (Mp,— L) are 
nonequivalent, due to anisotropic pseudodipolar inter- 
actions (described by the constant Jpd). 

In the absence of a field, states Al,2 and Bl,2 have 
the same energy and are equivalent. Stress and magnetic 
fields remove the degeneracy of these states, as can be 
seen from the phase diagram shown in Fig. 3. Solid 
lines in Fig. 3 represent the stability ranges of different 
states'* calculated by the direct minimization of potential 
(9). To calculate the stress values, we used an estimation 
[19] for the spontaneous strain uq — ^Mq/c44 « 10^^. 

Lines 2 and 3 correspond to the step-like (spin-flop) 
transition Bl,2— > Al (line 2) and A2 (line 3) accompa- 
nied by the 90° rotation of the magnetic vectors L and 
Mp. Within the regions that lie below line 1 and above 
line 2 (and below line 3 and above line 4) , potential (9) 
has only two minima that correspond to states Al and 
A2 with opposite orientations of the magnetic vectors. 
Such a "FM" state can be induced solely by a magnetic 
field. The effect of a stress reveals itself in the appear- 
ance of the two- or three-state regions (above line 1 (3) 
and below line 2 (4)), where equilibrium states Al and 
B (and A2 and B) have the perpendicular orientation of 
the magnetic vectors. It is also worth to note the quite 
general result mentioned in Sec. 2: the application of a 
stress can either enlarge or diminish (depending on the 
sign of cr) the stability regions of homogeneous states 



in-plane anisotropy X 10 



7.8 X 10- 



'* Since neither field nor stress in the accepted geometry discern 
between Bl and B2 states, we consider both states as B. 
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Fig. 4. (Color online) a — H phase diagram of inhomogeneous 
states for multiferroic Sr2Cu304Cl2. Left and right panels cor- 
respond to different field scales (below and above -ffdm)- Solid 
lines are calculated from the condition of ^j = 0, the composi- 
tion of equilibrium DS is represented schematically in each region. 
Double-head and single-head arrows show orientations of the AFM 
and FM vectors, respectively. ^''^•' = for simplicity 



and, in that way, invert the order of phase transitions 
between different states. 

Thus, the combination of a stress and a field opens 
a way to ehminate any state^ and, as a result, to pro- 
duce single-domain states Al (above lines 1 and 2) or 
A2 (below lines 3 and 4) ^^nd any combination of two- 
and three-domain states. 

In the next sections, we consider the behavior of dif- 
ferent multidomain states which appear in the finite-size 
sample under the application of a stress and a magnetic 
field. 

5. Stability Diagram of a Multidomain State 

In a finite-size sample, the equilibrium distribution of 
the vectors L and Mp is, in principle, inhomogeneous 
due to boundary conditions on the sample surface. We 
consider the multidomain states that consist of homo- 
geneously ordered regions (domains) separated by thin 
domain walls. No inhomogeneous distribution of L and 
Mp within the walls is taken into account. However, we 
assume that the domain walls can move under pondero- 
motive, field-induced forces, so that the DS can reach 
an equilibrium, after the field is applied. The DS is 



® States Bl and B2 could be discriminated by a small misalign- 
ment of the magnetic field from [110] direction, as it was done, 
e.g., in Ref. [17]. 



thus unambiguously described by a set of angle variables 
{9j, (fij} (j = Al, A2, Bl, B2) and fractions S^j of domains 
of the j'-type. Obviously, '^^j=l- The number of free 
variables in the energy potential (9) depends upon the 
types of movable domain walls present in the sample. 
The equilibrium DS in the presence of an external field 
is then found from the condition for the minimum of </> 
with respect to independent variables. 

In this section, we consider the case where all four 
types of domains can freely grow or diminish in size (say, 
in a virgin sample that initially contains domains of all 
types). If a stress a and a magnetic field H are rather 
small (region Al + A2 + B in Fig. 4), both fields are 
screened by an appropriate domain configuration, and 
the effective fields inside the sample are zero. Formally, 
this means that ^j are the only free variables. The equi- 
librium values of the magnetic variables in this case are 
given by Eq. (10), and the domain fractions depend on 
the magnetic field and the stress as follows: 



Ca1,A2 — -7 



1 _ ^(0) ± 
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Hdci 
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where we have introduced the notations 



i?dm = mpiV.^'^A/o, i?dcs = SVJo^Vd^, 



2A^' 



dcs 



CTO = 



C44 



\Ml 



i 



(0) ^ ^''2an 
^des " 



(11) 



(12) 



(13) 



Physically, i/dm is the FM demagnetizing field calcu- 
lated as if an AFM subsystem is absent. In an analogous 
way, -ffdos and ctq could be treated as the destressing field 
in the absence of FM ordering expressed in the mag- 
netic and stress equivalents, respectively. The shape- 
dependent value ^(°) represents the disbalance between 
type A and type B domain fractions in the absence of a 
field. For the sample under consideration (see Table 2), 
^(0) = 0.22 [19]. 

As seen from Eq. (12), the value of magnetic destress- 
ing field is enhanced due to the exchange interactions 
(constant Jq). On the contrary, the demagnetizing field 
is weakened due to a small FM moment (mp ^ 1). 
So, the demagnetizing field in the crystal under con- 
sideration is much smaller than the destressing ones. 
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^dm ^ ^des (see Table 2). However, for the chosen ge- 
ometry of the sample, the demagnetizing and destressing 
effects are of the same order of value, as will be shown 
below. 

The analysis of Eqs. (11) shows that the type of 
favorable domain depends upon the sign^ of ct/ctq. If 
a I oq > 0, domains of B type are unfavorable and disap- 
pear (^B — 0), when (see line 1 in Fig. 4)) 



equilibrium fractions 



ctb = ctq 



1+e 



(0) _ 



( H 



\Hdc 



(14) 



A negative stress (t/ctq < at i? = can remove domains 
of A2 type; respectively, ^ai = at (line 2 in Fig. 4)) 



0'A2 = CTq 



^(0) _ 1 + 



2H 

Ha,„ 



H 

-ffdns 



(15) 



Lines 1 and 2 in Fig. 4 bound the four-domain region 
and intersect at the point H = iJdm, a ~ ctq (l + f^°^). 
The two- domain state Al + A2 consists of FM domains 
and is stable in the region bounded by lines 1 and 4 (the 
latter corresponds to H — ifdm)- 

Three-domain state A1+B1+B2 (or Al+B, as Bl and 
B2 are equivalent) is equilibrium in a wide region below 
hues 2 (Eq. (15)) and 1 (Eq. (14)). The corresponding 
domain structure can be treated as an AFM one, because 
there are domains with different orientations of the L 
vector that compete with one another. It is worth to note 
that the negative stress a/ao < makes the domains 
of B-type favorable. The magnetization vector Mp of 
such domains is oriented perpendicularly to an external 
magnetic field H, as can be seen from the expression for 



Tab 


1 e 


2. Parameters used in numerical 


simula- 


tions. 


The source of data (exp 


3rimental or calculated) 


is speciflec 


in the last column. 


Destressing parameters 


corresponc 


to T = 120 K. 






Parameter 


Meaning 


Value 


Ref. 


a X b 


X c 


Sample size 


7 X 2 X 0.5 mm^ 


[17] 


Mo 




Sublatt. magnetization 


27.4 Gs 


[16] 


m-p 




Mp/Mo 


7x 10-* 


[16] 


Mp 




Satur. magnetization 


7 X 10^'^ emu/g 


[17[ 


^(0) 




Shape-induced bias 


0.22 


[19] 


MO 




Spent, strain 


10-6 


[19[ 


^dm 




Demagnetization field 


0.3 Oe 


Eq.m 


^dos 




Destressing const. 


7 mOe 


[19[ 


^dcs 




Destressing field 


2.1 kOe 


[19] 













^ The sign of tro oc A is material-dependent and, in the case under 
consideration, is unknown. 



Ui 


1 

~ 2 




/ ctq 


5 


&1 


^^B, 


1 
^ 4 


1 1 £(0) "^F^ 


^ H y 




V ^dcs / 



(16) 



In other words, the action of a stress may compensate 
the action of a field. A negative stress can even remove 
Al-type domains (with Mf||H) from the sample. This 
takes place at the critical line (line 3 in Fig. 4) : 



CTAl = CTq 



j(o) 



mpH 

4_^dos 



H 



(17) 



Two more things are also worth noting. First, the lines 
within the discussed phase diagram of the multidomain 
state correspond to the reversible (second-order type) 
transitions realized through the motion of domain walls. 
In contrast, the irreversible step-like switching between 
domains (related to the disappearance or nucleation of 
a certain state) can occur while crossing the stability 
lines within the phase diagram of the homogeneous state 
(Fig. 3). 

Second, the scale of the reversible phase diagram 
(Fig. 4) is proportional to a small shape-dependent de- 
magnetizing field ifdm — 0.3 Oe and a destressing stress 
ctq oc 40 kPa. The scale of the irreversible phase dia- 
gram (Fig. 3) is defined by the intrinsic spin-flop field 
Hs-i oc 1 kOe [19] and the characteristic stress scale 
oc 1 MPa and is at least two orders of magnitude larger. 
This scale difference is due to the different activation 
barriers and, hence, different kinetics for the nucleation 
and the domain wall motion. 

In the next section, we consider the behavior of the 
DS consisting of two types of domains with different ori- 
entations of L vectors. 

6. Competition of Two Domains: Nontrivial 
Field Dependence of Strain 

Let us consider a sample which was preliminary mon- 
odomainized to state Al by the excursion into the re- 
gion of high fields (above lines 1 and 2 in Fig. 3). 
Suppose that the magnetic field H is then diminished 
at a fixed stress value, a. Which kind of the DS 
can be expected in this case? It was already shown 
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// 
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U 



-200 -150 -100 -50 50 100 150 200 

H.Oe 

Fig. 5. (Color online) Field dependence of the macroscopic strain 
m(H). The sample contains domains of two types, Al and B. 
Lines 1, 2 and 3 correspond to cr = 0.4,-0.4, and —2 kPa, respec- 
tively. Inset shows the same field dependence for a single domain 
(B-type) sample. Stress and magnetic field are applied as shown 
schematically in the right bottom angle 

in [19] that, a,t H > Q and cr = 0, metastable ho- 
mogeneous state B has lower energy than homoge- 
neous state A2. The same relation is true for a neg- 
ative stress, Act < 0, as the energy of state B in this 
case is diminished, and also takes place in a rather 
wide range of positive stress values (between lines 2 
and 3 in Fig. 3). So, the nucleation of domains 
of B type is much more probable than the nucle- 
ation of A2-type domains. If, in addition, there is a 
slight misalignment between the magnetic field H and 
a crystal axis [110] that removes the degeneracy be- 
tween the Bl and B2 states, the DS of the sample 
is represented by the domains of only two types, Al 
and Bl. 

To illustrate the peculiarities of macroscopic proper- 
ties of such two-domain state, we calculated the de- 
pendence of the macroscopic magnetization, M™'*'^™, 
and the elongation, u = (n-un — Tru/2) (where u is 
the strain tensor), on the magnetic field H at a fixed 
external stress. Equilibrium values of the variables 
{dj,(pj,£,j}, j = Al,Bl, in this case were calculated by 
the numerical minimization of potential (9), by using 
the data from Table 2 with limitations ^a2 — ^B2 = 0. 
The corresponding dependences are given in Figs. 5 
and 6. 

The particular feature of the dependence u(H) in 
Fig. 5 is its linearity which is due to the presence of 
Al domain with small FM moment. Really, in the 
single domain Al, the spontaneous strain is constant. 



0.8 ^ 


H^maoro 
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Fig. 6. (Color online) Field dependence of the macroscopic magne- 
tization Af ™'"^'^° (H) . The sample contains domains of two types, 
Al and B. Lines 1, 2 and 3 correspond to cr = 0.4, —0.4, and 
—2 kPa, respectively 



In the single domain B m(H) shows the quadratic, 
symmetric (with respect H — > — H) behavior (in- 
set in Fig. 5) that agrees well with the symmetry 
of the strain tensor with respect to the time inver- 
sion (and the corresponding inversion of H). How- 
ever, in the two-domain case, it is the domain frac- 
tion ^ that depends almost linearly on iJ, thus lead- 
ing to a nontrivial and unusual field dependence of the 
strain. 

It can be also seen from Fig. 5 that the external stress 
induces a bias (shift) of the u(H) dependence, which 
is due to the nonequivalence of domains Al and B in 
the presence of a stress. The same bias is induced into 
magnetization curves (see Fig. 6). The value of bias is 
controlled by the stress. In thin films, such a bias can be 
introduced into the sample either by an external stress 
or by a properly chosen misfit between the film and the 
substrate. 

It is important to emphasize the role of a de- 
stressing energy in the field-induced behavior of the 
strain and the magnetization. In the case un- 
der consideration, the demagnetizing (stray) energy 
is too small to ensure the formation of metastable 
B-type domains. Really, the demagnetization field 
ffdm = 0.3 Oe is much smaller than the characteris- 
tic value of magnetic field, at which the appearance 
of B-type domains is favorable, iJcr oc N'^^^/mp ex 
100 Oe ]19] (see also Figs. 5 and 6). So, the mag- 
netic field-induced switching between Al-type and Bl- 
type domains is possible only due to destressing ef- 
fects. 
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7. Conclusions 

We have considered the behavior of the AFM+FM muhi- 
ferroic under the action of two external fields of different 
nature. Our calculations show that such a combination 
of fields makes it possible to separate the field influence 
on the different coexisting order parameters. This, in 
turn, opens a way to control the DS and macroscopic 
properties and to produce the states with any desirable 
types of domains. 

The simultaneous application of two fields can also 
increase the susceptibility of a sample to one of the 
fields, like it is the case with the magnetic suscep- 
tibility of a strain. The range of the field values, 
in which the system is sensible to external fields, 
can be controlled by the appropriate choice of the 
sample shape and the corresponding shape-induced 
"de-" fields (demagnetizing, depolarizing, destressing, 
etc.). 

For the present work, we have chosen the multifer- 
roic with a rather simple, tetragonal structure. How- 
ever, it would be interesting to generalize the pro- 
posed approach to multiferroics with ferroelectric and 
magnetic order parameters that show simultaneously 
electro- and magnetostriction (like PZT-PFW systems 
[20]) or/and have incommensurate AFM structure (like 
BiFeOs) that can be realized in more than four types of 
domains. We believe that three types of possible con- 
trol fields (electric, magnetic, and stress) in the first 
case and a large number of possible equilibrium do- 
mains in the second case can give rise to some new 
peculiarities in the macroscopic behavior of multifer- 
roics. 
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TEOPETHHHHH OnHC 

<l)A30B0i /IIArPAMH MyjIbTHOEPOIKIB 
CKIHHEHHHX P03MIPIB Y 3MIHHHX MEXAHIHHE 
HAnPy>KEHH51-MArHITHE nOJlE: KOHKyPEHIlM 
MI>K <DEPO- TA AHTHOEPOMAFHITHMH flOMEHAMH 



BHH, B 5IKHX cniBiCHyiOTb flBa THHH BnopjiflKyBaHHH, MO:>KHa Kepy- 
BaTH 3a flonoMoroio nojiiB pisHo'i npiipo^H. B .z^aHiii po6oTi npo- 
anajiisoBaHO noBeflinKy MyjibTHCJ^epoiKiB 3 aHTHCJ^epo- (AOM) Ta 
(J^epoMarniTHMM (OM) BnopsiflKyBaHHHM nifl flieio soBHiniHBoro 
MarniTHoro nojiH Ta MexaHi^Horo Hanpy>KeHH5i. KoM6iHaLi,i5i li,hx 

flBOX nOJliB /],03BOJI5ie OTpHMaTH MaKpOCKOniMHi CTaHH 3 pi3HOIO 
flOMeHHOK) CTpyKTypOiO. OTpHMaHHH TaKHM CnOCo6oM .Z],BOflOMeH- 

HHH CTaH fleMOHCTpye HCTHnoBy j^Jin AOM jiiHiiiHy 3ajie>KHicTB 
MaKpocKoniiHoi fle(J)opMari;ii Bifl 30BHimHboro MarniTHoro nojiH. 
Majie, ajie HenyjiBOBe MexaniMHe Hanpy>KeHH5i npiiBO.z],HTi> TaKo:>K 
flo 3cyBy KpHBOi 3ajie>KHOCTi MaKpocKoni'^HOi HaMarHinyBaHOCTi 
Bifl MarHiTHoro nojia (e4)eKTHBHOMy "nifliviarHiHyBaHHio") . 



O.B. FoMOHau, G. F. KopnieHKO, B.M. Jlonmee 

P e 3 lo M e 

MaKpocKoniMHHMH BjiacTHBOCTHMH Myjii>TH4)epoiKiB , to6to peMO- 
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